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Abstract 

The internal motion in a molecule, in which isomerization processes oc- 
cur, is characterized by two essentially different modes of motion - oscillatory 
and rotational. The quantum equation of motion which describes an isomer- 
ization process is reduced to the Mathieu-Hill equation. As is known, this 
equation is able to describe both modes. In the paper, it is shown that the 
chaotic region of the energy spectrum characterizing an isomerization process 
corresponds to the region where two modes of motion are merged together. 

Keywords: 

Isomeric reactions, dynamical stochasticity, nonlinear dynamics, quantum 
chaos 

PACS: 82.30. Qt, 05.45.-a, 05.45.Mt 



1. Introduction 

Dynamic stochasticity [l], [sf occurs in the region of parameter values, 
where topologically different trajectories adjoin each other. Trajectories near 
the boundary are highly sensitive to small disturbances. For small distur- 
bances these trajectories take a very complicated shape, which is one of the 
manifestations of stochasticity. 
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An elementary example of such systems is a pendulum that may perform 
two kinds of motion - oscillatory and rotational. This explains in the main a 
growing interest in physical processes in which motion reduces to pendulum 
type motions. 

A quantum analogue of dynamic stochasticity - quantum chaos 0, Isf 
occurs in the same region of parameter values where dynamic stochasticity 
does. In that case, the quantum states of two different symmetries adjoin 
each other. A highly complicated process takes place in a neighborhood of 
the boundary - the energy spectrum corresponding to the symmetry of one 
kind transforms to the energy spectrum of the other kind, which, in turn, 
demands the removal of degeneration corresponding to the symmetry of one 
kind and the appearance of degeneration corresponding to the symmetry of 
second kind. In the energy spectrum graph this process is represented by 
the set of branch points and merged energy levels. In passing through these 
points, the system may lose information on the initial conditions and pass 
from the pure state to the mixed one. The mixed state thus formed will be 
the state expressing quantum chaos. From that moment of time, the state 
of the system is not any longer described by a wave function and it becomes 
necessary to describe it by a density matrix 

Freque ntly , the objects of quantum chaos realization are polyatomic molecules 



id . I 111 . Il2[. For example, quantum chaos may take place in a polyatomic 
molecule, one of the generalized coordinates of which performs torsional os- 
cillations which, with an increase of the amplitude, transform to rotational 
motion [12]. Quantum chaos is a characteristic feature of such a transitional 
region. 

From the standpoint of quantum chaos it is of special interest to inves- 
tigate the so-called floppy molecules which have a capacity to isomeriza- 
tion j^, 10|. An example of such molecules is lithium cyanide (LiNC). This 
molecule contains a firm fragment C = N with a triple bond and a relatively 
light lithium atom (Li) which is able to perform angular oscillations relative 
to the molecule axis. When the oscillation amplitude attains the threshold 
value, the lithium atom detaches itself from the nitrogen atom and passes 
over to the carbon atom (Fig. la). At higher energies, it again returns to the 
nitrogen atom and there begins the alternating jump-over process, i.e., the 
lithium atom begins to rotate about the fragment C = N. 

The energy spectrum of lithium atom motion was studied by the numer- 
ical methods in jo], [lo[. It was shown that the energy spectrum consists 
of equidistant and nonequidistant regions. The energy levels arranged in 
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a nonequidistant manner correspond to lithium atom states at which the 
oscillatory mode of motion transforms to the rotational one. 

The aim of the present paper is to study the energy spectrum of lithium 
atom motion relative to the firm fragment C = N. As different from the 
papers j^, which are based mainly on the numerical methods, our ap- 
proach to the problem is analytical: using numerical data we approximate 
the potential energy of the lithium atom and study the energy spectrum by 
the analysis of the Schrodinger equation. 

2. An Isomerization Potential 

The lithium atom motion relative to the firm fragment C = N can be 
described by means of two coordinates R and 6 only (Fig. lb). 6* = 0° and 6 = 
180° correspond to the equilibrium states of Li-CN and NC-Li, respectively. 
R is the distance between the lithium atom and the centre of mass of the 
fragment C = N. 




Figure 1: The izomerization reaction in lithium cyanide: a) lithium atom motion relative 
to the firm fragment; b) the coordinates R and describing the lithium atom. 

The potential energy V of lithium atom motion relative to the firm frag- 
ment C = N was studied in |9| (Fig. 2). It was shown that the potential 
energy has two different minima for ^ = O'^ and 6 = 180°. The energy dis- 
tance between these minima is F(180°) — V{0) ~ 6.21 kcal/mol, while the 
distance from the minimum at = 0° to the maximum is 3.42 kcal/mol. 
Totally, the isomerization energy, i.e. the distance from F(180°) to Vmax is 
^ 9.63 kcal/mol (Fig.2). 

It can be easily verified that the potential energy of isomerization is ap- 
proximated by the formula 

V{9) = C + ai cos 9 — a2Cos'^9, 
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Figure 2: The plot of the potential energy V{6) obtained by a numerical method (9|. 
The angle 9 is laid off along the abscissa axis, and the potential energy in atom units 
(1 am.u. = 627.52 kcal/mol, 1 kcal/mol = 7 ■ 10~'^^J) is laid off along the ordinate axis. 



where ai < la^- Indeed, the points of a minimum of the potential energy 
V^(^) are Q = ttu, while the points of a maximum 6^- = ± arccos ^ + 2?™. 
If «! — )■ 2^2 5 then 6-^ ^ 6- ^ 2Tm. This means that minima vanish at 
points 27rn, while maxima adjacent to them merge. Therefore ai < 2a2 is 
the condition of existence of two pits with different depth values. 

We begin the reading of the potential energy V{9) from the largest min- 
imum y[(2n + l)7r]. Then C = ai + 0^2, and for defining the numeri- 
cal values of ai and 02 we obtain the equations V{6) = 2ai, V^tc) = 0, 



v{e.. 



± 



{2a2+ai) which for the above-mentioned numerical 
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values we obtain a = 3.245 kcal/mol and ^2 = 6.243 kcal/mol. 

The plot of the approximated potential energy is shown in Fig.3. It should 
be noted that the existence of potential minima of different heights is typical 
of isomerization processes where the firm fragment is nonsymmetric 13| . 

A Hamiltonian function describing the lithium atom motion relative to 
the firm fragment C = N has the form 



(1) 



where / = MiE? + M2r^ is the moment of inertia, and Pe is the generalized 
pulse corresponding to the cyclic variable 6. 
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V(9) 




Figure 3: The potential energy of the hthium atom representing the isomerization process. 
The graph is constructed for the potential energy V{6) = a2 + ai + ai cos 6 — a2 cos^ 6, 
when ai < 2a2- 



111 111 

TT^ + ' TT^ + ' (2) 

Ml rriLi rriN + mc M2 uin mc 

V{9) = a2 + Oil + ai cos 9 — a2 cos^ 9. (3) 

iTT'Li, and mc are the masses of lithium, nitrogen and carbon atoms. The 
expression (1) does not take into account the radial motion of the lithium 
atom, which is assumed to be at some average distance from the fragment 
C = N. 



3. The Classical Consideration. The Phase Portrait 

Let us first consider the lithium atom motion relative to the firm fragment 
C = in classical terms. For the sake of simplicity, we will study this motion 
for energy E > 2a. As will be seen later, these energy values include the 
oscillation and rotation modes which we are interested in. 

To obtain the motion equation, we write the energy integral 

+ CKi + q;2 + COS 9 — a2 cos^ 9 — E. (4) 

Using this integral of motion we can obtain a period of the lithium atom 
rotation about the fragment C = N 
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where 



de 



27 - U{0) 



7± 
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1 ± 



«2 + 1 



2a2 ) J 



E 
a.2 
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F{a,k) 



ai 
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0L2 



\/l — k"^ sin^ (p 



is a eliptic function of first order. 



When E +E. 



(2a2+ai) 
4a2 



in the hmit from above we have a — > 0, 



7± — > |, r — > 1, and the rotation period 
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2^2 V2' 
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In 



— >■ oo. 



Therefore as E ^ +Ec the rotation period logarithmically tends to infinity. 
Analogosously, we can calculate the oscillatory motion period T_. These 
calculations are not given here. We only note the oscillatory motion period, 
too, diverges logarithmically when E — >■ —Ec, i.e., if £■ — >■ —Ec, then T_ — >■ 
-Tc oo. 

As a rule, the existence of a logarithmically diverging period indicates the 
presence of two different modes of motion and the existence of a boundary 
trajectory - separatrix between them. This means that the separatrix appears 
to be obtained at energies E — —Ec- 

Using the expression (4) we can construct the phase portrait on the plane 
{9, 9). As seen from Fig.4, the phase portrait consists of closed (elliptic) tra- 
jectories, which correspond to the oscillatory process, and wave trajectories 
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describing rotational motion. The alternation of elliptic trajectories of two 

different dimensions on the ^^-axis is caused by the presence of two differ- 
ent minima of potential energy. Thus we see that lithium atoms may have 
two types of motion - oscillatory when E < ^^"4^°^^ and rotatory when 

E > ■ As seen from the figure these two types of motion are divided 

by the separatrix. 




~^ Zj I I I I I I I I I I i_J I I I I I I I I i_J I lZ a 

-10 -5 s 10 

Figure 4: The phase portrait corresponding to an isomerization process, constructed by 
means of the equation (4). 



4. Quantum-Mechanical Consideration of Motion. The Mathieu- 
Hill- Schrodinger Equation. 

Let us now consider the motion of the lithium atom relative to the frag- 
ment C = N from the standpoint of quantum mechanics. We replace the 
Hamiltonian function (1) by the Hamilton operator 

H^^ + V{0), (5) 

where p = —ifi^ is a pulse operator. Then the Schrodinger equation takes 
the form 

- i^'oQ^ + l*^! + 0:2 + ai cos 9 - a'i cos^ 6')'0 = £■0- (6) 
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We make the replacement 6 ^ 26 and introduce the notation 



a2 

«2 + y 



qi = ^ai, 



Q2 



2h? 



a2. 



Then, by simple transformations, from the equation (6) we obtain 



+ [E- 2X{qi cos 26 + cos 40)]^ = 0. 



(7) 



For the sake of simplicity, we restrict our consideration to the investigation 
of the dependence of energy levels on one (and not on two) parameter and to 
this end we introduce the parameter A into (7). This relationship is studied 
assuming qi and q2 to be given values. The energy spectrum of the initial 
problem is obtained at A = 0.5. The Schrodinger equation takes the form 
known as the Hill equation which in turn implies the Mathieu equation as 
qi — )■ 0. Using the known methods ij], solutions of this equation are sought 
for in the form of Fourier series 



Ce2n{z;X;qi,q2) 
Ce2n+i{z] X]qi,q2) 
Se2n+i{z; X;qi,q2) 



J2A?;\X)cos2rz, 



r=0 



Y,Atl\'\X)cos{2r + l)z. 

r=0 

oo 

r=0 

oo 



£';t'^(A)sin(2r + l)^, 



Se2n+2{z;X;qi,q2) = ^ ^^"^^^(A) sin(2r + 2)^;. 



(8) 



r=0 



Substituting these values into the equation (7) and equating the coefficients 
of the same harmonic, we obtain recurrent relations for the coefficients A2r, 
A2r+i, B2r, -B2r+i, (their upper indexes are omitted). For example, for the 
first series of (8) these relations have the form: 



r 
r 
r 
r 
r 



> 



aAf) — XqiA2 — Ag2^4 = 0, 

1 (a - 4)^2 - 2XqiAo - XqiA^ - Xq2A2 - Aga^e = 0, 

2 (a - 16)^4 - XqiA2 - AgiAg - 2Ag2^o - Ag2^8 = 0, 

3 (a - 36)^6 - XqiAi - XqiAs - Ag2^2 - Ag2^io = 0, 

4 (a - 4r^)y42r - XqiA2{r-l) - XqiA2{r+l) - Ag2^2(r-2) - 
-Ag2^2(r+2) = 0. 



(9) 
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This is an infinite system of equations tliat establislies a relation between tlie 
coefficients A2r- In order that this system would have a nontrivial solution 
it is necessary that Det{Aij) = 0, where 

Aj = [a- 4(z - 1)^] 6ij - Xqi6i+i-j - Xq2Si+2j - Mi^i-j+i - Ag2^i;j+2 - 

-XqiSi-2Sj-l - Xq2Si;3Sj-i - Ag2^j;2^j;2, i = 1 . . . oo, j = 1 . . . oo.(lO) 

are elements of the matrix composed of the coefficients of the system of 
equations (9) and written by means of the Kronecker symbol. Analogously, 
we obtain recurrent relations for the coefficients A2r+i, -B2r and -B2r+i, matrix 
determinants and characteristic equations, the numerical solution of which 
gives the energy spectrum. 



5. A Quantitative Estimate of an Energy Spectrum 

To obtain a quantitative estimate of the energy spectrum, we have to 
calculate first the numerical values of the parameters qi and q2 contained in 
the Mathieu-Hill equation (7). 

Applying the obtained values of ai and a2, taking into consideration the 
carbon, nitrogen and lithium mass values {mc = 19.92 • 10~'^'^ kg, mi^ = 
23.24 ■ 10"^"^%, rriLi = 11.62 ■ 10~'^''^kg) and also the values of the linear 
parameters r and R from ^ (r = 2.186ao = 11.57 ■ 10~^^m, R = 4.05ao = 
21.43 ■ lO^^^m,), we obtain, by means of the equations (2),gi = j^ai ~ 9300 
and q2 = —%Ci2 ~ —8942. Using the obtained values of qi and q2 and 
the numerical solutions of the characteristic equation, we can construct the 
energy spectrum as a function of Mathieu characteristics E{X) (Fig. 5). 

As seen from Fig. 5a, for low energies we have levels arranged equidis- 
tantly, which corresponds to oscillations of the lithium atom near the nitro- 
gen molecule, i.e., to oscillatory motion in a deep {6 = 0°) pit (Fig.3). When 
the lithium atom energy is e > 6 kcal/mol, the alternating transition from 
one pit to the other and vice versa takes place. In Fig. 5 to this transition 
there correspond levels arranged nonequidistantly. The spectrum obtained 



by us qualitatively coincides with the picture obtained in [10|. In both cases 



nonequidistant levels appear in the same intervals of energy. 



6. Conclusion 

Finding of energy spectrum of states which corresponds to the process of 
isomerization comes to the finding of proper values of Mathieu-Hill equation. 
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Figure 5: The spectrum of lithium atom motion about the fragment C = N . Spectrum 
given in the Fig. 5a is obtained from the Fig. 5b by means of projection on the energetic 
axis of the values of energy at A = 0.5. a) we see the nonequidistance of energy levels 
of the lithium atom at energies £ > 6 kcal/mol. In this region, the transition from the 
oscillatory mode to the rotational one takes place, b) the energy spectrum depending on 
the parameter A. 



The energy spectrum of lithium atom motion includes a region e < 
4kcal/mol, where levels are arranged equidistantly. These levels correspond 
to small oscillations of the lithium atom near isomeric states. Levels ar- 
ranged nonequidistantly begin for relatively high energies, when the energy 
level reaches half of the barrier height and exceeds e > Qkcal/mol. 

It should be noted that the system considered in the paper is integrable. 
Therefore the generality of nonequidistant (chaotic) levels is not the man- 
ifestation of quantum chaos though it may appear when these states are 
disturbed by periodic force as illustrated by the example of a quantum pen- 
dulum 0, i- 
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